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The largest eigenvalue of the reduced density matrix for quantum chains is shown to have a simple 
physical interpretation and power-law behaviour in critical systems. This is verified numerically for 
XXZ spin chains. 



In a recent paper Eisert and Cramer [l| investigated the following problem from quantum information 
theory: Given an initial state of a total system, what maximally entangled state between two subsystems 
can be obtained from it by local transformations in both parts, together with classical communication ? 
The general answer is that the spectrum of the reduced density matrix for the new state must majorize 
that of the initial one [2|. If the final state is maximally entangled, it is the superposition of M product 
states which all have the same weight, and the reduced density matrix p then has M non-zero eigenvalues 
with magnitude 1/M. To reach this state, the largest density-matrix eigenvalue w\ in the initial state must 
therefore be smaller than 1/M. This lead the authors to study this eigenvalue and the associated quantity 
S\ = — Inn?! for free-fermion chains. For critical systems, they found a logarithmic behaviour S\ ~ InL if 
the length L of the smaller subsystem diverges. This is the same behaviour as for the entanglement entropy 
S = — tr(plnp). Moreover, the factor of proportionality for the XX chain was 1/6 which is just half the 
value one finds for S. In the following we show that this can be understood very easily and generalize the 
result to arbitrary conformally invariant models. 

Consider an infinitely long quantum chain in its ground state and an interval of L consecutive sites. The 
reduced density matrix can be written in the form 

p=iexp(- J ff) (1) 

where Z = tx(exp(—H)) such that p is normalized to one, tr(p) = 1. Then the largest eigenvalue of p is 
given by 

wi = -exp(-Eo) (2) 

where Eg denotes the smallest eigenvalue of H. Therefore one has 

S 1 =\xiZ + E (3) 

Now consider 

In tr(>") = In tr(cxp (-nH)) -n\nZ (4) 
in the limit of large n. Then only the smallest eigenvalue of H contributes to the first term and one has 

In tr(p n ) = -nE -n\nZ (5) 
Therefore the quantity Si is given by the limit 



S 1 = - lim -lntr(p") 

n — *-oo fi 



(G) 
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By comparison, the entanglement entropy can be expressed as 

S = \nZ+<H> (7) 

and follows from values n near one via 

S = -^tr{p n )\ n=1 (8) 

Thus both S and Si are determined by tr(p"), or equivalently tr(exp (— nH)). The latter quantity 
has a simple meaning in terms of the classical two-dimensional system associated with the quantum 
chain 0, For this system, exp(— H) is the partition function for a plane with a cut of length L 
along which the variables are held fixed. Thus tr(exp (— nH)) gives the partition function of a system 
on a n-sheeted Riemann plane which has the topology of a double winding staircase. Therefore Si is 
the difference between the dimensionless free energies of one section of the Riemann surface and of a 
normal plane. The entropy S, on the other hand, is obtained from an almost planar system of conical 
shape, and its logarithmic behaviour can be traced back to the critical free energy of such systems 0, IB- ■ 

To determine 5*1 for a critical system, one can now use the conformal result Q 

c, 1, 
L — in ) 

tx{p n ) = b n {-) 6 V n J (9) 
a 

where c is the central charge in the conformal classification and a denotes the radius of small circles around 
the branch points which one has to exclude in the continuum model. The constant b n is not determined 
from the conformal considerations but independent of L. Inserting @ into © one obtains immediately 

s l = c -M-) + ki (io) 

o a 

where a non-universal constant k\ has been added. The logarithmic term is exactly one-half of the one 
appearing in S. For the XX chain, c = 1 and one recovers the result found in In the same way one can 
obtain Si for the other situations which have been studied. If the subsystem forms the end of a half-infinite 
chain, the prefactor of the logarithm is reduced by 2 and if one considers a chain of length L which is divided 
into two sections of length I and L — I, the result is Q 

r IT irl 

5 1 = -ln[ — sin(-)]+fc 1 (11) 
12 ira L 

Although the conformal results give the general answer to the problem, it is still instructive to look at the 
free-fermion case once again. Then the operator H in JQ) is also fermionic and has the diagonal form 

H = ^e k c\c k (12) 

k 

where the single- particle eigenvalues Ek follow from the correlation functions in the considered state H,EJ. 
For a half-filled chain of fermions, corresponding to an XX model, these quantities have been studied before 
|l0llll| and appear in pairs (e, — e) if L is even. Then one can express S in the form 



£ fc >0 



ln[l + exp(-e fc )] ' 



exp (e fe ) + 1 



(13) 
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whereas Si is given only by the first sum. For odd L, an additional eigenvalue Ek = appears and gives a 
contribution In 2 to both S and Si. This is valid for arbitrary L. In the limit L — * oo, however, the situation 
simplifies as the spectrum of the Sk then becomes linear and dense. For even L it is given by 

£fe = ;rrr (2fc + i); * = 0,1,2,... (14) 

2 IxiL 

One can then go over to integrals in l|13|) and finds that the two contributions to S are exactly the same 
(and proportional to hiL). This gives directly S± = S/2 for the leading logarithmic terms. The argument 
can be extended to non half-filled systems, where the spectrum of the Ek is shifted. 



In the free-fermion case one can also make use of the results for the Ek obtained for non-critical XY spin 
chains and for the Ising model in a transverse field 0, Il2l Il3| to calculate Si . The same holds for the XXZ 
model for A > 1. As the entanglement entropy, Si is then finite and diverges only as one approaches the 
critical point. Using the formulae in [l2T |. it can be expressed in closed form. For example, a transverse Ising 
chain divided into two halves gives in the disordered region 



Si 



where k measures the coupling, k' 2 = 1 — 
criticality (k — > 1) this has the expansion 



1 

24 



ln(- 



16 



)-7T 



I(k') 



(15) 



k 2 k' 2 ' " I{k) 

k 2 and I(k) is the complete elliptic integral of the first kind. Near 



Si 



1 

24 



ln( 



1 -k' 



ln( 



1 



(16) 



The logarithms can be expressed in terms of the correlation length, since £ ~ 1/(1 — k). This expansion 
differs from that for S, where the second term does not appear. For the XY chain one obtains similar results. 



Returning to wi itself, the results imply that it varies as a power in l/L for critical systems. Its vanishing 
for L — > 00 is connected with the fact that the spectrum of p, as well-known from DMRG calculations, 
flattens more and more as L becomes larger. The power-law behaviour can also be seen in numerical 
calculations. This is shown in Fig.l where results for XXZ spin chains with anisotropies A between —1 and 
+1 are presented. They were obtained from DMRG calculations on open chains of odd total length up to 
801 sites, keeping m=500 states but using no sweeps. The division was into subsystems of L and L — 1 
sites with L odd. Thus one is dealing with the situation described by the formula (fill) . The odd length was 
chosen because then the convergence is better. As one can see, the curves in the double-logarithmic plot 
become rather straight for large L and differ by additive constants. The resulting values of c are shown in 
the right part of the figure as function of 1/lnL. They clearly tend to the value 1, but the convergence is 
relatively slow and dependent on the value of A, as well as on the chosen expression for the slope of Si. 
A second order extrapolation in 1/lnL gives c- values which deviate from 1 by at most 5%. By contrast, 
the values for c determined from the entanglement entropy S lie much closer to 1 and converge faster. One 
could associate this with the different behaviour of the two quantities near criticality mentioned above. 

Summing up, we have shown that the largest density-matrix eigenvalue for a quantum chain can be 
interpreted as a difference of two free energies. Correspondingly, it shows universal features if the system is 
critical. For the problem cited at the beginning, this means that the rank M of the maximally entangled 
state grows unlimited with the subsystem size (although the power is relatively small) and the central 
charge c appears again in the field of quantum information. 

Note added: The conformal result for Si was given independently also by Orus et al.0] with an additional 
subleading term. 
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FIG. 1: Left: Largest eigenvalue of p for XXZ chains with half-length L. Right: Central charge c deter- 
mined from the slopes in the left figure multiplied by 12. The curves correspond to the anisotropies A = 
-0.9; -0.5; -0.1; 0; 0.1; 0.5; 0.9, from top to bottom. 
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